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ABSTRACT 

Error  estimates  for  Gaussian  quadrature  are  given  in  terms  of  the  number 
of  quadrature  points  and  smoothness  properties  of  the  function  whose  integral 
is  being  approximated.  An  intermediate  step  involves  a  weighted-L1  polynomial 
approximation  problem  which  is  treated  in  a  more  general  context  than  that 
specifically  required  to  estimate  the  Gaussian  quad  ature  error.  <  — • 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  gives  error  estimates  for  Gaussian  quadrature,  one  of  the  most 
widely  used  methods  for  approximating  integrals.  It  also  studies  a  related 
approximation  problem  that  has  more  general  potential  applications  to  the 
study  of  numerical  approximations  .  A  particular  anticipated  use  of  the 
results  on  Gaussian  quadrature  is  for  error  estimates  of  the  discrete 
ordinates  method  for  the  transport  equation,  by  the  second  author  in 
collaboration  with  Professor  J.  Pitkaranta  of  the  University  of  Helsinki. 

The  transport  equation  appears  as  a  model  equation  in  nuclear  engineering, 
astrophysics  and  climatology,  to  name  only  a  few  instances. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC ,  and  not  with  the  authors  of  this  report. 


ERROR  ESTIMATES  FOR  GAUSSIAN  QUADRATURE 
AND  WEIGHTED-L1  POLYNOMIAL  APPROXIMATION 

Ronald  A.  DeVore  and  L.  Ridgway  Scott 

1  .  Introduction 

""-e  purpose  of  this  note  is  to  establish  error  estimates  for  Gaussian 
quadrature  that  reflect  the  fact  that  the  accuracy  is  not  degraded  by  certain 
singularities  of  the  function  to  be  integrated  at  the  boundaries  of  the  inter¬ 
val  of  integration.  While  it  is  well  known  that  polynomial  approximation  can 
achieve  greater  accuracy  at  the  boundary  than  in  the  interior  (cf.  Timan  [3], 
p.  262),  error  estimates  for  Gaussian  quadrature  of  the  explicit  type  given 
here  have  apparently  not  been  presented  before .  The  basic  thrust  of  these 
estimates  is  to  show  that  the  error  in  N-point  Gaussian  quadrature  for  approx¬ 
imately  integrating  f  on  1-1,1]  is  bounded  by 

(1.1)  C  N~s  J1  |f<s)(x)|(1  -  x2)3/2dx 

8  -1 

for  all  integers  s  <  2N  such  that  above  integral  makes  sense.  Here  Cg  is 
a  constant  independent  of  N  and  f . 

Our  technique  of  proof  is  in  two  parts.  Firstly,  we  analyze  a  Peano 
kernel  for  the  quadrature  error.  This  has  the  effect  of  establishing  the 

bound  (1.1)  for  s  *  1  and  reducing  further  estimates  to  a  weighted-L1 
approximation  problem  for  polynomials.  This  part  of  the  analysis  is  presented 

in  section  2.  Secondly,  we  consider  weighted-L1  approximation  by  poly¬ 
nomials  and  prove  estimates  for  the  error.  This  is  done  in  section  3. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  Nos.  MCS-81-01661  and  MCS-81-01631 . 


Finally,  the  estimates  of  sections  2  and  3  are  oombined  in  section  4  to  give 
error  bounds  for  Guassian  quadrature  of  the  type  (1.1)  for  all  s  *  1.  The 
calculations  in  section  2  are  done  specifically  for  Gaussian  quadrature, 
although  they  clearly  extend  to  more  general  quadrature  approximations  based 
on  orthogonal  polynomials  for  weighted  integrals. 

The  weighted-L1  approximation  estimates  show  that,  for  a  rather  general 
class  of  weights  w  >  0  and  sufficiently  smooth  f  ,  there  is  a  polynomial 
P  of  degree  at  most  N  so  that 

(1.2)  P  I  f  ( x)  -  P(x)|w(x)dx  <  C  n“8  pH  (9>(x)|(l-x)8/2w(x)dx 
-1  8  -1 

for  any  s  such  that  N+1  >  s  >  1.  In  the  case  w  =  1  and  s  »  1,  (1.2) 

has  been  established  by  N.  X.  Ky  [1] . 
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2 .  Error  estimates  for  Gaussian  quadrature 

Consider  Gaussian  quadrature  approximation  of  the  form 
1  N 

(2.1)  /  f (x)dx  -  l  U)  f(x  )  »:  I  (f) 

-1  j-1  3  3 

<  (x^)  are  the  zeros  of  the  Legendre  polynomials  and  {tuj  are  che  integrals 
of  the  associated  Lagrange  interpolation  polynomials,  cf .  G .  Szego  [2]).  The 
ordering  -1  <  x^  <  * • •  <  x^  <  1  will  be  assumed,  and  we  introduce  Xg  := 

-1  and  Xfl+i  :*  1.  We  wish  to  establish  estimates  of  the  error 

(2.2)  e  ( f )  :=■  f1  f(x)dx  -  I  (f) 

N  ^  N 

in  terms  of  N  and  properties  of  f.  For  example,  the  Peano  kernel  theorem 
allows  us  to  write 

e  (f)  «  J1  K(t)f ' (t)dt, 

-1 

at  least  for  smooth  f,  where  K(t)  =  e^fHj.),  |t|  <  1,  (cf.  (3.3)  below) 
and  Hfc  is  the  Heaviside  function 


(2.3) 


JO  x  <  t 

\  1  x  >  t. 


It  follows  that 

K(t)  -1-t-I  to  -  I  to  -  t  -  1 . 

X  ,  >  t  3  X  <  t  3 

3  j 


The  Chebyshev-Markov-Stieltjes  inequality  (cf.  G.  Szego  [2],  p.  50)  implies 
that 


1 


+  x  <  1  (i> 

J  i*1 


‘  1  +  Vi'  3 


I,1 


1  ,N. 


Therefore,  for  j  =  1,***,N, 


Xj_1  -  Xj  <  K(x^-)  «  0  <  K(x^+)  <  xJ+1  -  x 


Let  K  denote  the  continuous,  piecewise-linear  function  (with  respect  the 
knots  x^ ,  *  *  *  ,x^j)  which  interpolates  0  at  ±1  and  the  value 


s=  maxfx^  -  x^.x^  -  x^}  at  x^ , 
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j  =  !,•••, N 


Since  K  ia  itself  piecewise  linear, 

|K(x) |  <  K(x)  for  all  x  6  [-1,1]. 

To  bound  the  values  dj»  recall  (cf.  G.  Szego  [2],  p.  122)  that 

(2.4)  Xj  -  -cos  8^  where  (2N+1)0j/TT  e  [2j-1,2j],  j  -  1, •••,«. 


Thus 


*j  ’ 


j-1 


008 


j-i 


0080  , 


0 

«  /  ^  sin  0  d0 

9j-1 

(2.5) 

<  (8 j— 0 j  1 )  max  Csin  ® j_1 »  ain 


<  (3»/2N)  max tain  ain  ® j ^ * 

Since  the  function  0  ♦  sin  0/3  is  decreasing  on  [0,*],  (2.4)  implies  that 

V4sin  6,  ,<  sin  9  <  4  sin  6  ,  2  <  j  <  N. 

j-1  j  j-1 


Thus 


d^  <  6w  sin  Q^Af 


6* 


■C 


oos20j  /  N 


6* 


/N. 


/  2 

Since  the  function  x  ♦  *1-x  is  concave,  this  implies  that 

(2.6)  |K(x)  I  <  K(x)  <  6»/l-x2  /N  for  xe  [-1,1]. 

Thus  we  have  proved  that 

1  e  ( f )  I  <  6*N  1  f'  |f'  (x)  |/l-x2  dx. 

N  -1 
Since  (2.4)  implies  that 

1/N  <  sin  */(2N+1)  <  /l-x2 

for  x  6  [x  ,x  ]  and  N  >  2,  (2.6)  implies  further  that 

l  N 

(2.7)  |K(x) |  <  6*(1-x2) 

for  x  9  [XjtXg].  But  for  x  9  [-1,Xj]  1x^,1],  |K(x)  |  -  1  -  |x|.  Thus 


(2.7)  holds  for  all  x  6  1-1,1],  giving  the  following  estimate: 


These  estimates  hold  for  all  f  e  L1 ( [—1,1] )  whose  weak  derivative  is 

2 

integrable  with  respect  to  the  weight  1-x  as  can  be  seen  by  approxi- 

s 

mating  f  via  smooth  functions  (cf*  the  definition  of  the  space  Y  given 

w 

later  in  section  3).  Summarizing  the  above,  we  have  the  following: 

Theorem  1 »  Let  eN(f)  denote  the  error  in  N-point  Guassian  quadrature 
applied  to  f  e  L^([-1»1])  (see  (2.1-2)  for  definitions).  If  the  weak 
derivative ,  f ' ,  of  f  is  Integrable  with  respect  to  the  weight  1  -  x  , 
then 


If  f 1  is  integrable  with  respect  to 


|eN<f)l  *  N~  ^  dx* 

Estimates  involving  higher  derivatives  of  f  can  also  be  derived  by 

estimating  Peano  kernels.  For  example,  one  may  write,  for  any  1  <  k  <  2N, 

eN(f)  =  /'  Kk(x)f{lt)(x)dx 

^  ( V-2 ) 

(K1  =  K  in  the  previous  notation).  is  a  C  '  piecewise  k-th 

degree  polynomial  (with  knots  {x^})  satisfying  0  =  K^_(±1)  =*  •••  = 

I k— 1 )  f  k  1  k  (k— 1 1 

(±1)  and  =  (-1)  between  knots.  Moreover,  ' (x^+) 

(k-1 ) 

-  K/  (x  -)  =  o)  for  j  =  1, • • *,N.  Using  these  facts  together  with  a 
k  3  j 

special  oscillation  property  of  Kj,  it  can  be  shown  that 
|K  (x)|  <  C  min{( 1-x2) /N2 , ( 1-x2) 2). 

However,  it  becomes  increasingly  difficult  to  estimate  the  higher  degree 
kernels  K^,  k  >  3.  Thus  the  'cl lowing  approach  proves  more  fruitful. 

Observe  that  eN( f )  =  eN(f  -  P)  for  any  P  @  *2N-1'  where  denotes 

the  set  of  polynomials  of  degree  not  exceeding  r.  Therefore,  Theorem  1 
implies  that 

I e  ( f ) I  <  inf  f1  If'-Pl(x)  w(x)dx, 
peP  -1 

2N-2 


J  2 

/1-x  ,  then 
2 
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Am  r 


(2.8) 


1 

I 

where  w  is  the  weight  function 

w(x)  »  6*  min{/l-x2/N,  1-x2}. 

Thus,  estimates  of  «N( f )  are  reduced  to  a  weighted-L1  approximation  problem 

for  polynomials.  Such  problems  will  be  considered  in  the  next  section.  The 

results  of  section  3  on  approximation  will  be  combined  with  the  results  of  the 

present  section  in  section  4  to  give  higher  order  estimates  of  the  form 

(2.9)  |e  (f) |  <  C  !'  |f(8)(x)fmin{(/l-x2/N)“,  (1-x2)*)dx 

N  8  -1 

for  arbitrary  positive  integers  s  <  2N. 


\ 

\ 


i 

I 
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3.  Weighted-!^  approximation 


In  this  section,  we  shall  prove  a  weighted-L1  polynomial  approximation 
result  of  the  form 

(3.1)  inf  /’  |u-P| (x)w(x)dx  <  C  N-8/1  |u (S> (x) |w(x) ( 1-x2) S/2dx, 

pePs  -1  8,W  -1 

where  PN  denotes  the  set  of  polynomials  of  degree  not  exceeding  N,  s  is  a 
positive  integer,  C  is  a  constant  independent  of  N  and  u  ,  and  w  is 

S  *  W 

a  positive,  integrable  weight-function  to  be  discussed  in  more  detail  subse¬ 
quently.  Rather  than  state  our  abstract  conditions  on  w  and  u  initially, 
we  shall  develop  them  in  the  course  of  deriving  (3.1).  However,  suffice  it  to 
say  that  (3.1)  will  be  proved  for  a  class  of  weights  including  the  Jacobi 
weights 

(3.2)  w(x)  -  (1+x)a(1-x>e»  a,B  >  -1. 

As  a  first  step,  recall  the  Heaviside  function  Ht(x)  defined  in 

(2.3) .  If  u  is  sufficiently  smooth,  we  may  write 

(3.3)  u(x)  *  P (x)  +  t  J1  H  (x)(x-t)S  1u^8^(t)dt 

'  V  3- I  J  I  ^ ^  t 

for  all  x  e  [-1,1],  where  P,  e  P_  . .  Let  X  e  P  be  an  arbitrary 

i  a  i  t  N— S^l 

family  depending,  say,  piecewise  continuously  on  t,  and  define 

(3.4)  P(x)  *  P^(x)  +  -  -  J1  Xt(x)(x-t)S  1u^S,(t)dt. 

Then  P  e  PN,  and  Holder's  inequality  and  Pubini's  theorem  imply  that 

J1  |u-P| (x)w(x) dx 
-1 

(3.5) 

<  TTrTTl  r  t/1  |Xfc-Htl (x) |x-ttS  1w(x)dx} |u(s) (t) |dt. 

Thus  to  prove  (3.1),  it  suffices  to  construct  X  in  such  a  way  that 

(3.6)  J1  |X  -H  |  (x)  |x-t|  8  V(x)dx  <  C  N  *(  1-t2 )  8/,2w(t)  . 

-1  t 

For  t  near  ±1,  this  is  relatively  easy  to  <3o  under  the  following 
assumption  t 
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(A1 )  There  is  a  constant  A1  euch  that 

i)  w(x)dx  <  A^{1-t)w(t)  for  0  <  t  <  1 

and  ii)  /*  w(x)dx  <  A  ( 1  +t) v»{ t)  for  -1  <  t  <  0. 

-1 

a 

Example  1 .  Assumption  (A1)  holds  for  the  Jaoobi  weight  w(x)  *  (1+x)  (1-x) 
provided  a,0  >  -1 . 

Proof .  To  see  this,  it  suffices  to  verify,  say,  1) .  Then 

Pw(x)dx  =  /1(1+x)a(1-x)0dx 
t  t 

<  2m*xia’01  f\l-x)*dx 

t 

<  2max{o*o}(i_t)1+^/(i+3) 

<  [21<X,/(1+@)]  (1+t)a(l-t)1+6 
-  A^l-Owtt).  // 

Lemma .  Suppose  assumption  (A1)  holds.  For  1  “  0,1,  let  X1-  denote  the 
functions  given  by 

Ai(x)  -  i  for  all  x  e  1-1,1]  . 

Then  choosing  Xfc  «  X®  for  t  >  0  and  X  =*  X^  for  t  <  0  gives 


Proof .  Suppose  t  >  0.  From  (2.3)  and  (A1),  we  have 


r  I A  -H  | ( x) |x-t | S  ^(xjdx  =  I x-t | S  1w(x)dx 

-1  t 

<  (1-t)S  1  f^w(x) dx 

t 

<  a^i-oVd 

<  A^  ( 1-t2)Sw(t) . 


The  case  t  <  0  is  similar.  // 

Corollary.  Let  <  >  1  be  arbitrary,  and  suppose  (A1)  holds.  Then  for  any 

N  >  1  and  t  e  (-1,1)  such  that  /l  -t 2  <  ic/N,  the  choices  for  A  given 

in  the  Lemma  satisfy  the  estimate  (3.6)  with  C  =  A_kS. 

- s  I 

/  2 

In  view  of  this  corollary,  it  suffices  to  assume  that  /1-t  >  k/N, 

where  k  can  be  chosen  later  at  our  discretion.  We  shall  subsequently 

construct,  for  certain  values  of  t,  A  =  A  e  P  „  .  such  that 

t  t,r,N  N-2r+ 1 

(3.7)  IH^-A^lfx)  <  C^(/ l-t2/N) 2r  ^x-tl  2r  1  for  x  e  [-1,1], 

where  r  is  any  positive  integer  (to  be  chosen  later  depending  on  w  and 
s).  Furthermore,  A^  will  be  monotone,  nondecreasing,  with  A^(-1)  =  0  and 
A^(+1)  =  1.  Hence  also 

(3.8)  |H  -A  | (x)  <  1  for  all  x  e  f-1 , 1 j . 

Assuming  these  properties  of  A  for  some  t  for  the  moment,  we  proceed  to 
prove  (3.6)  for  such  t.  Let  S  =  /l-t  /N.  Then 
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(3.9) 


Z1  |H  -X  ||x-t|8"1w(x)dx 
-1  1  z 

<  /  |x-t|S  ^w(x)dx  +  C  /  |x-t|8  2lw(*<)dx  52r  1 , 

|x-t|<$  r|x-t|>6 

Ix|<1  |x|<1 

where  we  used  the  bounds  (3.8)  and  (3.7)  on  the  first  and  second  integrals, 
respectively.  Thus  the  estimate  (3.6)  follows  easily  from  the  following  two 
assumptions : 

(A2)  Let  s  be  a  positive  integer .  Then  there  exist  oonstants  A^  <  ®  and 

Y  >  2  such  that,  for  all  5  >  o  and  t  6  [-1,1]  satisfying  1-t2  >  Y$» 

/  w(x+t) |x|8_1dx  <  A,w(t)«8. 

lx |<«  2 

(A3)  There  exist  A^  <  “*»  Y  *  2  and  0  <  <  *  such  that,  for  all  6  >  0 

and  t  e  [-1,1]  satisfying  1-t2  >  Y$  and  all  k  >  kg, 

/  w(x+t)|x|  kdx  <  A_w(t)5  15+1  . 

I x | >6  3 

|t+x|<1 

In  applying  (A2)  and  (A3)  to  (3.9),  note  that,  for  6  <■  /l-t2/N ,  the  oondi- 
2  / - 2 

tion  1-t  >  Y^  is  equivalent  to  »1-t  *  Y/N  •  This  will  therefore  be 

satisfied  if  <  >  Y,  a  requirement  we  now  impose  on  «.  Choosing  r  such 
that  2r-s  »j  k  >  k^  thus  proves  (3.6)  for  t  such  that  (3.7-8)  hold. 

Before  proving  the  bounds  (3.7)  and  (3.8),  we  show  that  (A2)  and  (A3)  hold  for 
Jacobi  weights . 

a  0 

Example  2:  Assumption  (A2)  holds  Po c  t*ie  Jacobi  weights  w(x)  «*  (1+x)  (1-x) 
for  all  a, 6  e  R. 

Proof .  Take  Y  “  3.  By  a  change  of  variables, 

5  8  /  (w(x+t)/w(t) ) | X ! B  1dx 

I  X|<« 

r1- 

|x  |<0 

-  /  ( i  +  <7^r>y)a(1  - 

|y|<1  1+t 
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For  o  <  0,  one  has 

I  <  fM(  1-M  \)^y  kdy 
1 

<  /M(  1-m  2dy. 

i 

A  simple  calculation  shows  that,  for  M  >  1, 

/**(  1-M  V)Sy  2dy  <  CQ  for  all  8  >  -1, 

1  0 
completing  the  proof .  // 

Me  now  construct  X  satisfying  (3.7-8).  Mi  shall  do  so  by  choosing 

(3.10)  X  (x)  /*  8  (y)dy, 

c  -1 

where  5  e  P  approximates  the  Dirac  8-function.  In  particular,  8 

t  N-2r  t 

will  be  nonnegative  and  have  integral  on  t-1,1]  equal  to  one.  Thus  X^ 
satisfies  the  claimed  monotonicity  property  as  well  as  having  the  required 
values  at  ±1.  Hence  (3.8)  will  follow  automatically. 

To  oonstruct  8^,  first  note  that  it  suffices  to  take  8fc(x)  :*■  Vj  for 

/  2 

all  x  0  [-1,1]  for  N  <  4r  (recall  that  we  restrict  t  to  satisfy  »1-t 
>  Y/N  in  (3.7-8)).  For  N  >  4r,  let  n  be  the  greatest  integer  not 
exceeding  N/2r.  Mote  that  N/n  <  3r.  Let  Tn  denote  the  Chebyshev 
polynomial  of  degree  n,  and  let  (t^sl  <  i  <  n}  be  its  zeros: 

Tn(x)  oos(n(oos”1x)) 

t^  oos((i  -  V2  )w/n)  oos  »  i  «  1,***,n. 

We  shall  establish  (3.7-8)  for  t  6  (t^:1  <  1  <  n) .  Define 

(3.11)  8  (X)  c. [T  (x)/(x-t.)]2r, 

in  1 

where  c,  is  chosen  so  that  P  8  (x)dx  >1.  To  estimate  the  size  of  c^, 

-1  *i 

observe  that 

|T  (x)|  -  loos  n6|  >  V2  for  |n6  -  i*|  < 
n  j 

-12- 


;:V7^4* 


Thus,  writing  x  =  cos  8,  we  see  that  for  |n9  -  iw|  <  — 


|T  ( x)/  ( x-t  ^)  |  >  V2) 00 s8  -  cosS^I 


>  V2  I  cos  ( i  +  ■^)*/n  -  0038^1  1 

>  [  ( 5w/3n)max(sin(i  + -^)Vn,sin9^}]  ^ 


>  [ ( 5t/n) sin8^]  , 

because  sin9/9  is  decreasing  on  {0,^1  •  The  measure  of  the  set 
(x  =  oos8:|n9  -  iir  |  <  it/ 3 }  is 

cos ( ( i  -  j)ir/n)  -  oos((i  +  -j)Vn)  >  (2if/3n)min{sin (i  i  -jjn/n) 


>  (2n/9n)sin9i , 


Therefore , 


c.  =  /  [T  (x)/(x-t.)]  rdx 
i  n  1 


>  ((5Vn)sin9i)  “(2w/9n)sin9i  . 


(3.12)  cx  <  24(5irsin9i/n)2r’1  -  24(5*/l-t2/n) 2r_1  <  C. (/l-t2/N)2r-1 
To  complete  the  estimate  (3.7),  observe  that 


|H  -*  I  (x)  -  c 
i  i 


Therefore, 


I*  (y)<Jy  X  <  t. 
-1  i 


J1  $t  (y)dy  x  >  ^  . 

X  i 


|H  -X  |(x)  <  c  /2  y"2fdy 
i  i  I x-t ,  | 


<  ci|x-til  2r+1/(2r-1). 


Combined  with  (3.12),  this  proves  (3.7)  for  t  »  t^ .  Hence  (3.6)  is  now 
verified  for  t  ■  t^. 


I 

I 


For  the  general  case  t1+1  <  t  <  t*.  define  Afc  ... 

since  sinfl/6  i8  decreasing  on  [0,*], 

1  /  2 

3  lr' 1-t  /n  <  *  min{sin8i/Sine  }/n 


•  Mote  that, 
Ci 


*  /  +^sin®d8  »t  -  t 
6  i  i+1 


Therefore, 


<  3*  min{sin0^ ,ain91+^}/n  <  37f/l-t2/n. 
I  lRt“Xtl(x)  lx-t|S~1w(x)dx 


”  ^  l(*)tx“tla  ^w(x)dx 

*  /  !<x)|x-t|8  1w(x)dx 


+  Z1  |H  -A t.|(x)|x-t|S_1w(x)dx. 

-1  i  i 

This  first  term,  via  (A2).  ls  hounded  by 


/  1|x-t|8"‘,w(x)dx  <  A2(ti-t)8w(t) 


<  *2(ti“ti+1>8*tt) 


<  A,(3»/lV/n)*w(t) 


*  cr/,tyT-t2/N)’w(t) . 


(The  application  of  assumption  A2- is  valid  provided,  e.g., 

t^l  -  t1+1)  >  -  t1+1).  The  reader  may  easily  check  that  this 

holds  if  the  constant  <  mentioned  in  the  Hilary,  and  the  subsequent 
discussion,  is  chosen  sufficiently  large,  depending  on  r  and  y  .  This 
observation  applies  as  well  to  the  application  of  assumptions  A2  and  A3  in  the 


*  ■  v. 


•ext  set  of  inequalities.)  Since  (3.7-8)  hold  for  t^,  the  second  term  is 

/ - 2 

bounded  by  (5  :=  /l-t  /N) 


C  5 
r 


2r-1 


/ 


x-tP2ltrti+i|  1 


|x-t.'.  2r  +  1|x-t|S  ^(xjdx 


/ 


|x-t|<2|ti-t1+1 | 


|x-t|8  'w(x)dx 


<  cr  (2<5 ) 


2r-1 


|x-t|S  2rw(x)dx 


|x-t|>2|ti-ti+1 1 


+  2SA2lti-t1+1|Sw(t) 


<  Cr(2«)2r_1  A3(2|trti+1l)S_2r+1w(t) 


♦  2SA2lti_ti+i 


<  C  (/l-t2/N)Sw(t), 
r  /  8 


assuming  (A2)  and  (A3)  hold.  Therefore,  (3.6)  is  now  proved  for  all  t  @ 

[-1,1),  and  hence  the  estimate  (3.1)  is  established  for  smooth  u  and  for 

weights  w  satisfying  assumptions  (Al-3).  In  fact,  in  view  of  the  Lemma,  one 

/  2  s 

can  improve  (3.6),  and  hence  (3.1),  by  replacing  the  expression  (»1-t  /N) 

J  2  g  2  s 

by  min{(/l-t  /N)  ,(1-t  )  }.  To  extend  the  result  to  more  general  u, 
define,  for  positive  integers  s, 

(3.13)  *ul  “  If 1  |u*k) (x) |(1-x2)kw(x)dx, 

'  k-0  -1 

where  u^8^  is  interpreted  as  a  weak  derivative.  Define 

(3.14)  Y8  =*  (u  e  L,1  ( [-1 , 1 )  ) :  »ul  <  »}  . 

w  loc  w,s 
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g  »  # 

Then  Y  is  a  Banach  space  having  C  ([-1,1])  as  a  dense  subspace  .  Using 

W 

this  density,  we  arrive  at  the  following  theorem,  which  summarizes  our 
results. 

Theorem  2.  Let  w  be  a  positive,  integrable  function  on  [-1,1]  satisfying 

8 

assumptions  (Al-3),  where  s  in  (A2)  is  some  positive  integer.  Let  u  @  Y^ 
(see  (3.13-4)  for  the  definition).  Then  for  any  positive  integer  N  >  s-1. 


where  PN  denotes  polynomials  of  degree  not  exceeding  N  and  Cs  w  is  a 
constant  independent  of  N  and  u. 


Remark  1 .  If  u* s* (x) ( 1-x2) s^2  is  integrable  on  [-1,1],  the  above  estimate 
may  be  simplified  to  yield 

inf  J1  |u-P|  (x)w(x)dx  <  C  N  8  J1  |u*S)  (x)  I  ( 1-x2)  S/^2w(x)dx. 

Pe  P  -1  S'W  -1 

N 


Remark  2.  Assumptions  (A2-3)  imply,  in  particular,  that  w(t)  >  0  for  |t|  < 
1,  unless  w  =  0.  But  such  a  condition  is  necessary  for  an  estimate  such  as 
(3.1)  to  hold.  To  see  this,  suppose  that  w  is  continuous  at  t  and  w(t)  * 
0.  Choose  a  sequence  {u^}  of  smooth  functions  with  u^9*  positive,  sup¬ 
ported  in  (t  :  |t|  <  1/j}  and  J1  u*S*(t)dt  ■  1  (take  j  >  1/ ( 1— | t | )  for 

-1  3 

1  g«1 

simplicity).  Then  u^  converges  in  L  to  u(x)  Ht(x)(x-t)  /(s-1)l 
as  j  tends  to  infinity.  For  u  ■  u^,  the  right  hand  side  of(3.1)  will  tend 
to  zero.  However,  the  left  hand  side  certainly  will  not  do  so  since  u  4 
for  any  N. 


*Proof:  Define  ur(x)  u(rx),  0  <  r  <  1.  Then  as  r  -*  1,  uf  -►  u  in  Yw< 

Let  ur  be  obtained  from  i*r  by  mollifying:  ur  :»  ur  *  where 

supp{6  ^  »  [-e,e]  and  <5  has  integral  one  and  is'smooth.  Then  ur  E 

smooth  on  [-1,1]  for  e  £  1-r,  and  uf  -*  ur  in  Y®.  Now  let  1, 

£:=*  1-r  ■*  0  to  get  u„  r  •*  u  in  Y  .  ' 

r  i  t  w 


is 


Higher-order  error  estimates  for  Gaussian  quadrature 


In  this  section,  we  combine  the  results  of  previous  sections  to  prove  an 
estimate  of  the  form  (2.9).  From  (2.8),  we  know  that 
I e  ( f ) |  <C  inf  f1  |f'-P|(x)w(x)dx, 

/  2  2 

where  w(x)  :=*  rain{*1-x  /U,1-x  }.  It  is  easy  to  check  that,  if  a  collection 
of  weights  w^  satisfy  (A1),  (A2)  or  (A3),  then  so  (toes  the  weight  function 
minfw^}.  Thus  Theorem  2  applies  to  the  weight  w  since  it  is  a  minimum  of 
two  Jacobi  weights,  yielding 

|e  (f ) (  <  C  J1  |f (x) |min{ (/)-x2/N) 3  1 ,(1-x2)8  ^wfxjdx 


/ 1  |f*°* (x) |min((/l-x2/N)S,(1-x2)S}dx. 


we  summarize  this  final  result  as 

Theorem  3 .  Let  e^lf)  denote  the  error  in  N-polnt  Gaussian  quadrature 
approximation  to  the  Integral  of  f  oil  (-1,1]  (see  (2.1-2)  for 
definitions)  .  Suppose  that  (1-x^)sf (x)  (weak  derivative)  is  integrable 

g 

on  (-1,1] ,  i  .e .,  f  e  Y  where  s  is  any  Integer  such  that  1  <  s  <  2N . 
Then 


|e  (f)|  <  C  J1  |f <S) (x) |min{(/l-x2/N)8 , ( 1-x2 )S}dx, 
”  8  -1 

where  Cfl  in  independent  of  N  and  f . 


Remark  3 .  If  f (x) (1-x2)8^^  is  integrable  on  [-1,1],  the  above  estimate 
simplifies  to 

|e  (f)|  <  C  N  °  J1  |  f  (x)  |  { 1-x2)8,/2dx. 

N  8  -1 

This  Is  the  estimate  anticipated  in  (1.1). 


I 


t 

I 
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